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Executive  Summary 


Abstract 


The  work  performed  under  the  aforementioned  AFOSR  grant  includes  research  on  satelhte 
based  communication  networks.  It  consists  of  two  parts.  In  the  first  part  we  investigate  some  of 
the  arising  resource  allocation  and  handover  issues  in  non-geostationary  mobile  satellite  systems. 
We  propose  a  channel  classification  scheme,  in  which  available  carriers  are  partitioned  into  classes 
and  each  class  is  associated  with  a  range  of  propagation  delays  to  the  satellite.  The  key  idea  of  our 
approach  is  that  users  with  similar  propagation  delays  are  assigned  to  the  same  class  of  carriers.  The 
suggested  infirastructvue  results  in  better  resource  utilization  and  reduced  call  blocking  rate  and  can 
be  implemented  with  low  signaling  load.  In  the  second  part  we  consider  information  dissemination 
through  satelhte  broadcast  delivery  and  more  specifically  the  identification  and  treatment  of  the 
issues  that  arise  in  novel  multi-stage  broadcast  scheduling-caching  applications.  We  investigate 
the  design  of  the  satellite  schedule,  the  cache  management  policies  at  the  ground  stations  and 
the  terrestrial  schedules  that  arise  in  joint  fashion  in  the  context  of  our  problem.  Our  results 
demonstrate  that  such  a  two  stage  system  is  feasible  and  can  provide  more  efficient  data  deUvery 
compared  to  a  single  stage  system. 


1  Resource  allocation  and  handover  issues  in  non-geostationary 
mobile  satellite  networks 

1.1  Resource  allocation  in  MEO  mobile  satellite  neWorks 

In  the  near  future,  existing  terrestrial  radio  networks  are  envisioned  to  integrate  with  satellite  sys¬ 
tems  so  as  to  provide  global  coverage.  In  order  to  estabhsh  communication  for  non-hand-held  and 
hand-held  user  terminals,  the  radio  fink  design  must  allow  fioll-  and  half-duplex  operation  respec¬ 
tively,  where  the  latter  is  desirable  whenever  radiation  power  restrictions  are  imposed.  In  addition, 
due  to  user  mobility  and  volatility  of  the  wireless  medium,  sophisticated  resource  managpniPnt  is 
required,  so  as  to  enhance  system  capacity  and  reduce  the  likelihood  of  blocking.  An  inherent 
problem  of  the  satellite  link  is  propagation  delay,  which  can  lead  to  inefiicient  resoiuce  allocation. 
This  problem  becomes  more  important  in  medium-earth-orbit  (MEO)  satellite  systems,  which  are 
characterized  by  large  propagation  delays  and  large  intrarbeam  delay  variations. 

The  main  contribution  of  this  research  initiative  is  to  investigate  some  of  the  arising  resource 
allocation  and  handover  issues  in  non-geostationary  mobile  satellite  systems.  We  propose  a  chan¬ 
nel  classification  scheme,  in  which  available  carriers  are  partitioned  into  classes  and  each  class  is 
associated  with  a  range  of  propagation  delays  to  the  satellite.  The  key  idea  of  our  approach  is 
that  users  with  similar  propagation  delays  are  assigned  to  the  same  class  of  carriers.  The  sug- 
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gested  infrastructure  results  in  better  resource  utilization  and  reduced  call  blocking  rate  and  can 
be  implemented  with  low  signaling  load. 

1.2  System  model 

A  satellite  constellation  in  MEO  orbit  is  considered.  The  projection  of  satellite  on  the  earth  is  the 
sub-satellite  point.  The  M  are  classified  in  L  groups  Bi, . . . ,  Bl:  A  group  Bi  contains  all  equidistant 
beams  from  the  sub-satellite  point.  Beams  belonging  to  Bi  are  type-i  beams,  for  i  = 

Satellite  gateways  (GWs)  constitute  the  access  points  for  users.  We  use  the  terms  “call”,  “user” 
and  “user  terminal  (UT)”  to  refer  to  mobile  users.  An  important  parameter  that  characterizes 
the  UT-satellite  link  is  the  satellite  elevation  angle  6  with  respect  to  the  UT,  which  is  the  angle 
between  the  line  that  defines  UT  horizon  and  the  line  that  connects  the  satellite  to  the  UT. 

A  FDMA/TDMA  access  scheme  is  assiuned.  The  available  spectrum  is  divided  into  carrier 
frequencies.  A  TDMA  frame  in  a  carrier  has  6  time  slots,  each  of  duration  T^,  and  a  user  traffic 
burst  occupies  one  slot.  A  channel  is  perceived  as  a  distinct  carrier-slot  pair.  A  UT  transmits 
and  receives  once  in  a  firame  period.  We  assume  half-duplex  communication  with  non-overlapping 
transmission  and  reception  intervals  of  a  burst.  A  user  can  optionally  use  diversity  within  the  same 
carrier.  If  a  UT  transmits  and  receives  once  in  a  frame  and  there  is  a  guard  time  tg  «Ta  between 
transmission  and  reception  intervals,  3  slots  are  required  for  single-channel  half-duplex  operation 
and  6  slots  are  needed  in  case  of  diversity.  A  static  resomce  allocation  scheme  within  each  satellite 
is  adopted  and  a  set  of  carrier  frequencies  is  assigned  to  each  beam. 


1.3  The  delay  class  concept 

Transmission  (Tx)  and  reception  (Rx)  traffic  bursts  of  the  user  are  separated  by  a  time  offset,  which 
is  referred  to  as  Tx/Rx  burst  offset  The  relative  positions  of  transmitted  and  received  bursts  in 
a  window  of  3  slots  depend  on  Tx/Rx  burst  offset,  which  in  turn  depends  on  the  UT  propagation 
delay  Tp  to  the  satellite.  Users  in  different  positions  in  a  beam  require  different  Tx/Rx  burst  offsets 
to  maintain  non-overlapping  transmission  and  reception  intervals.  If  these  users  are  assigned  to 
the  same  carrier,  several  slots  will  remain  unoccupied.  In  figure  1  the  situation  with  a  single  or 
multiple  offset  values  in  a  carrier  are  illustrated.  For  a  single  Tx/Rx  biurst  offset,  efficient  call 
accommodation  is  achieved  by  packing  users  in  contiguous  time  slots.  However,  for  users  with 
different  delays,  the  system  resorts  to  multiple  burst  offsets  in  order  to  maintain  non-overlapping 
transmission  and  reception  intervals  and  thus  a  significant  amount  of  resources  is  wasted  (in  the 
figure,  slots  marked  with  “X”).  For  a  given  arrival  rate,  blocking  rate  is  increased,  since  fewer  calls 
are  accommodated  in  the  system. 

In  MEO  mobile  satellite  systems,  the  situation  above  arises  in  edge  beams  with  large  intra¬ 
beam  delay  variations  due  to  the  curvature  of  earth  surface.  In  order  to  circumvent  the  problem, 
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we  introduce  the  concept  of  delay  class.  For  a  nominal  delay  Tq,  we  define  a  delay  range  5Tp  > 
0  around  Tq.  Transmission  and  reception  bursts  of  UTs  with  propagation  delays  in  the  range 
[To  -  5Tp,TQ-\-5Tp]  are  accommodated  within  a  reference  window.  UTs  with  propagation  delays 
in  that  range  belong  to  the  same  delay  class  corresponding  to  To  and  should  be  allocated  to 
the  same  carrier  or  group  of  carriers.  Several  delay  classes  with  delays  To,*  must  be  defined  and 
class  i  corresponds  to  a  carrier  group  Gi.  Delays  of  users  assigned  to  carriers  of  Qi  must  satisfy 
To,i  5Tp  <Tp  <  To^i  +  STp.  Pictorially,  each  nominal  delay  To,i  corresponds  to  a  contour  (circle)  on 
the  earth  surface.  All  nominal  delays  represent  concentric  circles,  centered  at  the  sub-satellite  point 
Q  (figure  2).  A  contour  of  delay  To^i  comprises  points  on  earth  with  that  delay  to  the  satellite.  The 
contours  of  delay  To,i  ±  STp  form  a  zone  that  is  defined  as  delay  class  i.  Biursts  of  UTs  in  a  delay 
class  must  arrive  aligned  at  the  satellite  and  UT  [1].  Offset  values  for  a  specific  UT  are  derived 
by  comparing  UT  propagation  delay  Tp  and  nominal  delay  To^i  of  its  delay  class  and  the  offset  is 
proportional  to  Tp  —  To^i.  The  range  STp  is  derived  by  considering  upper  and  lower  bound  of  delays, 
such  that  accommodation  in  a  window  is  feasible.  It  can  be  shown  [2]  that  STp  =  Ts/4  —  tg/2. 

Delay  classes  have  certain  positions  with  respect  to  the  beam  pattern.  Figure  3  shows  a  projec¬ 
tion  of  one  quadrant  of  the  beam  pattern  on  a  two-dimensional  plane  and  the  relative  position  of 
delay  classes.  Beams  that  are  closer  to  the  footprint  edge  are  elliptical  and  elongated  and  therefore 
they  cover  a  wide  range  of  delays.  A  delay  class  serves  a  certain  set  of  beams,  namely  beams  of  the 
same  beam  type  Bj  due  to  circular  symmetry.  Beams  of  a  beam  type  may  be  covered  by  more  than 
one  delay  class.  For  example,  out-most  beams  axe  covered  by  three  delay  classes,  since  propagation 
delay  range  is  larger.  Intermediate  beams  may  be  served  by  one  or  two  delay  classes. 

Resources  consist  of  carriers  that  belong  to  a  pool  and  are  assigned  on  a  per  beam  basis  with 
dynamic  or  fixed  channel  allocation  (DCA,  FCA)  schemes.  Although  DCA  schemes  capture  satellite 
movement  and  traffic  variations,  we  consider  fixed  allocation  of  carriers  to  a  beam  for  illustrative 
reasons.  If  arising  calls  are  uniformly  distributed  in  the  area,  the  expected  amount  of  traffic  in 
out-most  elongated  beams  is  larger  and  call  blocking  is  increased.  The  delay  class  concept  can  be 
used  to  alleviate  this  problem,  by  assigning  more  delay  classes  and  therefore  more  carriers  to  these 
beams.  The  number  of  delay  classes  k  and  their  positions  depend  on  satellite  orbit  and  height, 
footprint  and  beam  size,  frame  structure,  traffic  burst  length,  and  guard  time.  Since  adjacent  delay 
classes  may  overlap,  the  maximum  residence  time  at  a  delay  class  overlap  area  is  another  important 
design  parameter  that  affects  delay  class  positions.  This  residence  time  represents  maximum  allowed 
tolerance  for  delay  class  transitions. 

Let  (f)  be  the  earth  central  angle  (figure  4).  Let  d  and  r  be  the  distance  and  delay  from  a  delay 
class  contour  to  satellite.  Assume  that  H  is  the  satellite  height  and  Re  is  the  earth  radius.  A 
“cup”  on  the  earth  is  the  area  defined  by  a  delay  class  contour  and  a  pole,  and  a  “zone”  is  defined 
by  delay  class  contours,  i  and  i  +  1.  Let  A  be  the  area  of  such  a  cup  or  zone.  The  minimum 
and  maximum  elevation  angles  0mm  and  6max  stem  from  visibility  conditions  arid  definition  of  the 
closest  delay  class  to  satellite  nadir.  Range  [Omim^max]  represents  the  portion  of  footprint  to  be 
covered  with  delay  classes.  Equivalently,  this  area  is  defined  by  central  angle  range 
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where  the  relation  between  angles  (j)  and  6  is 


=  (1) 

and  (pmin  =  <f>{Smax),  <l>max  =  4>{^min)-  The  following  algoritlun  finds  delay  class  positions  without 
considering  delay  class  overlap  regions. 


•  Step  1  :  Obtain  the  pair  {(j>min,0max)- 

•  Step  2  :  At  first  iteration,  determine  the  first  delay  class.  Compute  parameter  ha  = 
ReO-  —  cosOjnax)  (figure  4),  i.e.,  the  height  of  first  (close  to  nadir)  delay  class.  A  cup  is 
thus  defined  on  earth  surface  at  height  ho  and  find  the  area  of  this  cup,  A  =  2nREho.  Then, 
compute  distance  parameters, 

al  =  R%-{RE-hof  md4  =  al  +  {H  +  hof  =  H^  +  2iRE  +  H)hQ,  (2) 

where  do  is  the  distance  from  the  first  {j  =  0)  delay  class  to  satellite  and  find  delay  to  the 
satellite  tq  =  do/c. 

•  Step  3  :  At  jth.  iteration,  determine  the  {j  +  l)th  delay  class.  Select  constant  height  hj  = 
A/{27rKRE)-  Then,  compute  distance  parameters 

=  and  d?  =  o|+  +  =H^  +  2{RE  +  H)'j2hi.  (S) 

\  i=0  J  \  i=0  /  i=0 

and  find  delay  of  {j  +  l)th  delay  class  to  satellite,  rj  =  dj/c. 

•  Step  4  :  Repeat  this  procedure  for  j  =  1, 2 . . . ,  k  —  1. 


Heights  hj,  defining  {j  +  l)th  delay  class  are  selected,  so  that  the  area  of  each  zone,  defined  by 
heights  hj^i  and  hj,  is  equal  to  that  of  the  defined  cup.  The  algorithm  presupposed  a  given  number 
of  delay  classes  k.  If  that  is  not  given,  it  can  be  derived  as 


K  =  min  : 


cos 


-1 


(4) 


The  satellite  footprint  is  sequentially  covered  with  delay  classes  until  angle  Omin  is  reached. 


1.3.1  Delay  class  determination  and  resource  allocation 

We  study  the  problem  of  determining  the  serving  delay  class  for  a  call.  When  a  call  is  initiated, 
resources  are  assigned  to  it  after  a  call  request  message  that  contains  current  satellite  and  beam 
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identities  and  the  propagation  delay  to  satellite.  If  a  beam  or  a  delay  class  handover  occurs,  the  new 
delay  class  can  be  determined,  since  timing  to  the  current  satellite  is  maintained.  However,  after  a 
satellite  handover,  synchronization  is  lost,  since  satellite  synchronization  systems  are  independent. 
Although  a  satellite  handover  is  less  frequent  than  a  beam  handover,  it  can  certainly  occur  since 
satellite  footprints  move  fast  on  the  earth  surface.  A  satellite  handover  may  occur  during  the  call, 
if  a  call  is  either  of  long  duration  or  is  carried  by  an  edge  beam.  The  derivation  of  the  delay  class  of 
the  new  satellite  is  important  for  reliable  resource  assignment  and  the  delay  from  satellite-UT  delay 
needs  to  be  determined.  Two  methods  can  be  used  by  the  GW  to  determine  this  delay.  According 
to  method  1,  the  GW  can  retrieve  the  most  recent  estimate  of  UT  position  and  associate  it  with 
the  new  satellite  ephemeris  data  to  find  an  estimate  of  the  delay.  In  method  2,  the  GW  can  request 
an  explicit  measurement  report  from  the  UT.  The  first  method  is  faster  and  easier  to  implement, 
while  the  second  one  is  more  accurate  but  also  more  time-  and  bandwidth-consuming.  The  first 
method  should  be  given  priority  and  used  whenever  the  estimated  UT  position  is  accurate  enough 
to  provide  a  good  estimate.  A  UT  position  error  is  acceptable  if  it  does  not  lead  to  an  incorrect 
delay  class.  If  the  UT  position  estimate  is  not  accurate  enough,  the  system  should  resort  to  the 
second  method.  The  question  that  arises  is  when  to  use  each  method,  so  as  to  minimize  the  incurred 
signaling  load.  One  could  argue  that  Global  Positioning  System  (GPS)  would  facilitate  accurate 
computation  of  delay.  However,  in  the  case  of  a  satellite  handover,  delay  measurement  with  GPS 
would  require  message  exchahges  between  UT  and  several  satellites,  which  would  consume  power 
and  bandwidth.  In  addition,  delay  determination  methodology  must  maintain  some  compatibility 
with  second  generation  technologies  that  do  not  include  GPS  in  their  standards.  Hence,  the  delay 
can  be  determined  with  the  two  aforementioned  methods. 


1.3.2  UT  position  error  tolerance  region 

Each  UT  is  characterized  by  an  actual  location  on  earth  and  a  delay  to  a  satellite.  Only  estimates 
of  these  quantities  are  available  and  estimated  UT  position  is  referred  to  as  known  UT  position. 
Consider  a  beam  that  is  served  by  two  delay  classes  (figure  5).  Let  the  delay  classes  be  defined  by 
delays  ro,i  ±  6Tp,  for  i  =  1,2,  with  To^i  <  To, 2-  When  the  UT  is  located  in  regions  1  or  2,  it  is 
assigned  to  a  carrier  of  the  corresponding  delay  class  (1  or  2).  Region  3  corresponds  to  delay  range 
[To, 2  -  STpyTo^i  +  STp]  and  is  the  overlap  region  of  delay  classes.  Range  6Tp  determines  the  length 
of  overlap  region  and  is  usually  fixed  for  any  pair  of  delay  classes.  While  in  the  overlap  region,  a 
UT  can  be  served  by  carriers  of  delay  class  1  or  2.  The  best  quality  carrier  of  two  carrier  groups 
is  selected.  The  UT  may  belong  to  one  of  the  three  depicted  regions  in  a  beam.  However,  because 
of  delay  estimation  errors  it  may  seem  to  reside  in  a  different  region  from  the  actual  one.  As  is  ^ 
explained  in  [2],  an  incorrect  delay  class  assignment  occurs  if  the  difference  between  actual  and 
known  position  corresponds  to  delay  difference  that  is  greater  than  the  length  of  overlap  region. 
Thus,  the  length  of  delay  class  overlap  regions  in  the  new  beam  after  satellite  handover  must  be 
determined.  For  a  beam  with  n  delay  classes  and  n  —  1  overlap  regions,  finding  the  minimum 
length  overlap  region  and  converting  it  to  distance  corresponds  to  computing  the  worst  case  error 
in  position  determination  for  which  a  correct  delay  class  assignment  is  derived.  The  correct  delay 
class  must  be  assigned,  so  that  UT  is  assigned  to  an  appropriate  carrier.  The  following  methodology 
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is  proposed  to  solve  the  problem. 


•  Step  A  :  Divide  the  set  of  satellite  beams  into  subsets  Si,  62,  •  • . ,  so  that  each  beam  in 
Bj  has  Kj  delay  classes.  Beams  of  subset  Bi  are  type-i  beams  and  form  a  toroid. 

•  Step  B  :  For  each  delay  class  pair  (i,z  +  1)  of  beams  in  Bj,  find  lengths  of  overlap  regionss 
Xij^  as  Xij  =  To^i  -  To^(i+i)  -  2\STp\,  where  To,i  >  To^(i_|_i).  For  subset  Bj  with  kj  delay  classes, 
i  =  1, . . . ,  /Cj  “  1.  Then,  for  each  Bj,  select  Xj  =  min^  Xij,  to  account  for  minimum  length 
overlap  region  for  a  beam  with  many  delay  classes.  Then,  Xj  is  the  maximum  allowed  error 
between  actual  and  estimated  delay  for  beams  in  Bj,  for  correct  delay  class  assignment.  Thus, 
Xj  =  maxjt  \Tpj^  —  fpj  =  emaxj,  with  Tpj^,  Tp^  being  actual/estimated  delay  for  user  k  in  Bj. 

•  Step  C  :  Find  distance  tolerance  Adj  =  c  ♦  emaxj^  c  =  3  x  10®  m/sec  is  the  light  velocity. 

•  Step  D  :  For  each  subset  Bj  and  delay  class  i  =  1, . . .  «j,  compute  the  two  distance  extremes 
to  satellite,  =  c «  To,i  ±  Adj  and  corresponding  central  angles  with  the  law  of  cosine  [3]: 

2Re{Re  +  H)  ’ 


•  Step  E  :  Compute  radius  of  each  circular  tolerance  region  Lij  =  {Re/2) 


,  and  for 


the  beam  set  with  Kj  delay  classes,  select  Lj  =  minimi. Lij  to  account  for  worst-case  error. 


1*3.3  Delay  class  assignment 


We  describe  the  procedure  to  select  the  appropriate  method  (method  1  or  2)  and  assign  the  correct 
delay  class.  Delay  class  assignment  is  equivalent  to  estimation  of  UT  position.  The  region  in 
which  the  UT  resides  and  the  quality  of  position  estimate  determine  the  method  to  be  used.  We 
assume  that  delay  and  Doppler  frequency  offsets  are  available  and  the  estimated  UT  position  is 
derived  from  estimates  of  delay  and  Doppler  frequency  that  are  randomly  distributed  around  the 
actual  values,  according  to  a  Gaussian  distribution.  Let  UTact  and  UTknown  be  the  actual  and 
known  UT  positions  and  assume  that  the  beam  has  n  >  1  delay  classes.  The  known  UT  position 
corresponds  to  non-overlap  region  A;,  A:  =  1, ...  ,n,  or  to  overlap  region  £  between  delay  classes  £, 
£+1,  £  =  1, . . .  ,n  —  1.  In  the  first  case,  method  1  can  be  used  if  UT  position  error  is  less  than  the 
minimum  length  overlap  region.  In  the  second  case,  we  need  to  compute  ri  =  To,(^+i)  +  5Tp  —  fp 

and  r2  =  Tp—  (Tq,^  —  STp),  which  denote  distances  of  the  known  position  from  the  two  sides  of  the 
overlap  region.  The  procedure  to  make  the  delay  class  assignment  is  as  follows. 


•  Step  1  :  Derive  tolerances  Lj  for  each  beam  subset  Bj. 
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•  Step  2  :  For  every  UT  in  the  satellite  footprint,  execute  steps  2-5. 

•  Step  2  :  Compute  actual  propagation  delay  and  Doppler  frequency  oflFset.  Generate  delay 
and  Doppler  frequency  offsets,  based  on  Gaussian  distribution  and  derive  known  UT  position. 

•  Step  3  :  Compute  distance  between  estimated  and  actual  UT  position  and  determine  the 
beam  and  beam  type  Bm  of  known  UT  position. 

•  Step  4A  :  Known  delay  Tp  corresponds  to  a  UT  position  in  a  non-overlap  region  £.  Then,  if 
\UTact  —  UTknownl  <  Lmy  then  use  method  1  to  find  delay,  else  use  method  2. 

•  Step  4B:  Known  delay  fp  corresponds  to  UT  position  in  an  overlap  region  £,  Then,  if 
\UTact  —  UTknown\  <  min{ri,T2},  then  use  method  1  to  find  delay,  else  use  method  2. 

•  Step  5  :  Use  this  delay  value  to  assign  the  call  to  a  delay  class. 


If  parameters  ri,  r2  are  computed,  min{Ti,r2}  is  the  closest  distance  from  the  overlap  region 
boundary.  The  condition  \UTact  “  ^  min{ri,T2}  means  that  the  error  region  covers 

overlap  region  i  and  part  of  non-overlap  regions  ^  or  ^  +  1.  Then,  only  method  2  can  guarantee  a 
correct  delay  class  assignment.  Otherwise,  if  \UTact  —  UTknmvn\  <  niin{ri,r2},  the  error  region  lies 
entirely  in  overlap  region  £  and  method  1  can  be  used. 


1.4  Handover  issues  in  mobile  satellite  networks 

The  problem  of  resource  allocation  in  mobile  satellite  networks  is  inherently  related  to  handover  [4]. 
Each  time  a  handover  occurs,  a  resource  allocation  problem  needs  to  be  solved  and  new  resources 
need  to  be  assigned  to  the  user.  The  objective  is  to  accommodate  as  many  users  as  possible  in  the 
system  and  guarantee  acceptable  link  quality  for  each  user.  Then,  the  system  can  respond  better 
to  a  potential  sudden  load  increase  or  link  quality  deterioration  and  the  likelihood  of  blocking  a 
connection  request  is  minimized  [5].  Efficient  resource  allocation  and  handover  techniques  aim  at 
satisfying  quality  of  service  (QoS)  requirements  of  users.  The  contributions  of  this  research  are  (i) 
the  establishment  of  a  framework  for  reliable  handover  prediction  at  the  GW,  (ii)  the  consideration 
of  meaningful  criteria  for  selection  of  the  appropriate  satellite  and  beam,  whenever  the  user  can 
choose  among  several  satellites  and  beams  for  handover  and  (iii)  the  proposition  of  methodologies 
with  which  handover  signaling  load  is  reduced.  Several  alternatives  for  handover  prediction  and 
path  selection  are  proposed  and  evaluated  [6]. 


1.4.1  Beam  selection  for  measurements 


The  standard  procediue  of  beam  signal  level  monitoring  in  GSM  terrestrial  cellular  networks  is 
used.  Each  user  monitors  the  received  signal  strength  of  a  set  of  broadcast  common  channels 
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(BCCH),  where  each  BCCH  corresponds  to  an  adjacent  beam  of  a  satellite.  The  GW  or  satellite 
base  station  (SBS)  periodically  commands  the  user  to  measure  signals  of  all  visible  serving  and  non¬ 
serving  satellites  and  creates  a  list  of  beams  to  be  measured.  These  measurements  will  facilitate 
handover  decisions.  The  list  comprises  the  set  of  beams  of  different  visible  satellites  that  cover 
the  user  location,  as  well  as  a  set  of  approaching  beams  of  serving  satellites.  There  can  be  one  or 
two  serving  satellites,  depending  on  the  use  of  diversity.  The  beams  of  the  first  set  are  candidates 
for  satellite  handover,  while  the  beams  in  the  second  set  are  candidates  for  beam  handover.  The 
procedure  takes  place  at  call  origination  and  during  the  session  of  a  call  and  is  depicted  in  figure  6. 

The  problem  of  high  handover  signaling  load  becomes  important  in  non-geostationary  satellite 
systems,  which  are  characterized  by  small  beams  and  larger  number  of  handovers.  In  order  to 
efiiciently  utilize  scarce  bandwidth  resources,  signaling  load  should  be  reduced.  The  required  t.imp 
to  perform  measurements  and  ultimately  make  the  handover  decision  depends  on  the  size  of  the 
measurement  list,  which  increases  with  the  number  of  visible  satellites  and  Tvith  diversity.  Dming 
connection  setup,  the  objective  is  to  minimize  this  time,  so  that  transition  to  traffic  mode  occurs 
fast.  An  alternative  method  in  which  only  one  beam  from  each  visible  satellite  is  measiured  requires 
less  time  to  accompfish.  Although  such  a  method  is  desirable  during  connection  setup,  it  may 
result  in  a  suboptimaJ  path  during  traffic  mode.  Numerical  results  in  [6]  indicate  a  significant, 
improvement  in  measurement  time  and  associated  signaling  load  with  this  technique. 


1.4.2  Path  selection 


Path  selection  algorithm  is  an  inherent  component  of  resource  allocation  and  takes  place  after  beam 
selection  and  before  handover.  Each  entry  k  of  the  list  is  a  pair  of  satellite  and  beam  indices  (s^,  hk). 
The  list  is  modified  as  follows.  All  possible  combinations  {k,  1)  of  single  elements  are  created  and 
appended  to  the  list.  Entries  with  a  measurement  below  a  given  threshold  are  eliminated  since  they 
indicate  unreliable  connection.  The  list  is  ranked  based  on  a  preference  factor,  so  that  the  first 
entry  has  the  highest  preference  factor.  Each  entry  is  a  single  or  a  diversity  “path”  that  is  eligible 
for  resomce  allocation.  For  entry  k  the  preference  factor  P*  is  a  function  of  satellite  elevation 
angle  6  and  the  difference  of  the  received  signal  strength  from  threshold.  If  the  entry  is  a  diversity 
path,  the  preference  factor  depends  also  on  the  azimuth  separation  angle  (j)  between  satellites.  A 
path  with  high  satellite  elevation  angle  is  more  preferable  for  assignment,  since  the  likelihood  of 
call  blocking  due  to  some  obstacle  is  reduced.  Furthermore,  a  high  signal  strength  provides  better 
quality  and  reduces  the  probability  of  blocking  due  to  unacceptable  signal  level.  Finally,  a  big 
azimuth  separation  angle  is  more  preferable,  since  there  are  fewer  chances  that  both  paths  get 
corrupted  due  to  an  unpredictable  blockage.  Thus  for  the  more  general  case  of  a  diversity  path, 

where  A,  B  are  factors  that  captiue  the  relative  significance  of  these  criteria. 
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1.4.3  Handover  prediction  algorithms 

When  a  UT  enters  a  beam,  it  is  mapped  to  a  coordinate  system  [7],  until  a  time  interval  of  some 
length  (e.g.,  1  min)  is  found,  so  that  the  UT  resides  in  the  current  beam  at  the  beginning  of  the 
interval  and  in  a  different  beam  at  the  end  of  the  interval  (figure  7).  Handover  must  occur  sometime 
during  this  interval  and  predicted  handover  time  is  the  mid-point  of  the  interval.  Motivated  by 
the  fact  that  additional  handovers  lead  to  increased  signaling  load  and  delay,  we  present  a  new 
criterion  for  handover.  The  key  idea  is  that  the  user  should  be  forced  to  reside  in  a  cell  as  much 
as  possible,  so  that  handover  rates  are  minimized.  Upon  creation  of  the  list  of  candidate  paths  for 
transition,  no  preference  factor  computation  is  required.  Instead,  the  entry  with  the  beam  in  which 
the  mobile  is  predicted  to  stay  the  longest  is  selected  for  handover.  This  beam  may  or  may  not 
belong  to  the  serving  satellite.  This  criterion  involves  fewer  computations  and  measurements  than 
path  selection.  User  residence  time  for  each  beam  in  the  list  is  computed  by  standard  coordinate 
system  computations  without  the  need  for  elevation  angle  computations.  More  importantly,  no 
signal  strength  measurement  exchange  is  necessary. 

This  prediction  technique  assumes  stationary  UT.  If  the  UT  is  in  motion,  successive  handover 
predictions  may  lead  to  prediction  time  errors.  Beams  move  at  different  rates  relative  to  the  earth 
surface,  since  they  follow  different  paths  around  the  earth.  Edge  beams  move  at  lower  rates  than 
central  ones.  UT  mobility  has  significant  effect  on  beam  transitions  in  the  slow-moving  beams 
and  small  effect  on  fast-moving  ones.  Additional  results  on  the  effect  of  mobility  on  anticipated 
handover  times  can  be  found  in  [7]. 

For  the  satellite  handover,  the  serving  satellites  at  current  time  tc  must  be  identified  and  updated 
to  a  future  time  tf  with  satellite  ephemeris  data.  If  0c/t(^/)  ^  10®  or  Oawi'f^f)  ^  10®,  then  a 
handover  occurred  at  some  time  t*e(tc,t/).  We  apply  bisection  on  that  interval  and  stop  after  n* 
iterations,  when  —  tc,n*  <  W,  The  predicted  satellite  handover  time  is  ts  =  +  ic.n*)/'^ 

[6], 


1.5  Numerical  results 

A  MEO  mobile  satellite  system  with  the  specifications  of  ICO  satellite  system  was  simulated.  A 
satellite  has  163  beams  of  8  t3rpes.  A  UT  can  be  served  by  a  satellite  when  its  elevation  angle  9 
exceeds  Omin  =  10°,  which  is  typical  for  a  rural  environment.  Within  a  satellite  footprint,  delay 
class  positions  are  computed  and  concentric  rings  as  those  in  figure  2  are  defined.  In  order  to  serve 
beams  close  to  nadir  with  one  delay  class,  we  set  the  closest  delay  class  to  be  the  ring  at  0  =  60 
degrees.  By  applying  the  delay  class  determination  algorithm,  we  found  that  the  maximum  number 
of  delay  classes  that  ensures  beam  coverage  is  3  and  beams  are  covered  by  1,2  or  3  delay  classes. 
The  distance  tolerance  values  that  guarantee  correct  delay  class  assignment  for  beams  with  2  and 
3  delay  classes  were  found  to  be  142fcm  and  2Akm  respectively.  Position  determination  is  more 
sensitive  in  beams  with  3  delay  classes,  due  to  delay  class  proximity  and  small  length  of  overlap 
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regions.  UTs  are  covered  by  one  GW.  Satellite,  beam  and  delay  class  handovers  can  occur  for  UTs. 
In  satellite  handover,  a  UT  always  uses  one  satellite  as  long  as  it  remains  visible.  A  beam  handover 
occurs  when  a  UT  moves  in  the  coverage  area  of  another  beam  within  the  satellite.  A  delay  class 
handover  occurs  when  the  UT  moves  to  another  delay  class  within  a  beam.  Since  beams  or  delay 
classes  overlap,  a  beam  or  delay  class  handover  occurs  randomly  during  the  time  when  the  UT  is 
in  the  beam  or  delay  class  overlap  area.  We  simulated  an  one-hour  satellite  revolution.  Calls  in 
different  beams  arrive  in  Poisson  streams  of  rate  A  and  their  durations  are  exponentially  distributed 
with  mean  r  =  l//i  =  150sec,  which  is  typical  for  voice  transactions.  Traffic  intensity  is  measured 
in  Erlangs  (E),  E  =  Ar/60.  Although  handovers  are  affected  more  by  satellite  velocity,  a  simplistic 
UT  mobility  model  is  also  used  with  a  random  number  of  UTs  having  speed  uniformly  distributed 
in  [0,  72]  km/h  and  direction  in  [0, 27r]. 

The  carrier  grouping  method  leads  to  efficient  resource  management  and  reduced  blocking  rate. 
The  presented  mkhod  does  not  significantly  affect  blocking  rate  in  moderate  size  beams  with 
one  delay  class,  but  it  offers  a  clear  advantage  in  edge  beams  with  large  delays.  We  performed 
experiments  for  an  edge  beam  and  measure^  performance  for  randomly  generated  UTs  within  the 
beam  region,  assuming  that  carriers  are  allocated  to  the  beam  with  FCA.  In  the  scenario  with  no 
carrier  grouping,  UTs  were  allocated  randomly  in  carriers  irrespective  of  their  location  in  the  beam. 
Thus,  UTs  assigned  at  the  same  carrier  had  different  Tx/Rx  offsets  so  as  to  maintain  orthogonality. 
In  the  second  scenario,  the  beam  was  divided  in  3  delay  classes,  the  set  of  carriers  was  divided  in 
three  subsets  and  each  subset  was  mapped  to  a  delay  class.  Calls  belonging  to  a  delay  class  were 
assigned  to  a  carrier  of  that  class.  Each  time  a  carrier  of  a  delay  class  became  full,  UTs  were 
assigned  to  another  carrier  of  that  delay  class.  Overlap  area  calls  were  assigned  to  the  least  loaded 
carrier  of  those  in  the  delay  classes.  When  all  carriers  of  a  delay  class  were  full,  the  call  was  blocked. 
The  blocking  probability  is  the  ratio  of  blocked  calls  over  total  number  of  calls.  Experiments  with 
12  and  45  carriers  were  performed  and  results  are  illustrated  in  figure  8.  The  improvement  in 
performance  due  to  carrier  grouping  is  evident.  For  example,  in  the  case  of  12  carriers  and  a  traffic 
of  25  Erlangs,  blocking  probability  was  reduced  by  almost  50%.  For  higher  loads  or  for  a  larger 
resource  pool  of  45  carriers,  the  advantage  of  carrier  grouping  becomes  more  substantial. 

We  also  measured  the  resulting  satellite,  beam  and  delay  class  handover  rates  in  our  simulation. 
In  figure  9,  the  handover  rates  are  depicted  as  a  function  of  time,  for  a  square  of  15°  x  15°  within 
15°  longitude  and  latitude  distance  from  the  GW.  The  aggregate  handover  rates  for  all  traffic 
carried  by  the  GW  is  also  depicted.  A  first  observation  is  that  all  handover  rates  converge  to  a 
steady  state  after  a  transition  interval  of  about  2, 500  seconds.  In  steady  state,  about  30  delay  class 
transitions  per  minute  occur.  Delay  class  handover  rate  depends  on  whether  the  UT  is  located  in 
edge  beams  with  more  than  one  delay  classes.  The  ratio  of  satellite,  beam  and  delay  class  handover 
rates  in  steady  state  was  is  2  :  7  :  3.  For  locations  further  from  the  GW,  this  ratio  was  1  :  15  :  1. 
Furthermore,  delay  class  handover  rate  becomes  in  a  way  predictable,  since  it  occurs  only  within 
beams  with  more  than  one  delay  class. 

In  figure  10,  we  present  comparative  results  about  satellite  and  beam  handover  rates  by  using 
UT  position  or  maximum  beam  residence  criterion  in  a  region  with  moderate  load  (0.92  calls/sec). 
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With  the  second  criterion,  a  reduction  to  beam  handover  rate  up  to  85  —  90%  was  observed,  while 
for  heavier  traffic  load  this  reduction  reached  35%  ^  40%.  It  is  important  to  observe  the  increased 
satellite  handover  rate  for  some  time  periods.  However,  the  ejffect  of  increased  satellite  handover 
rate  is  low,  since  the  satellite  handover  rate  at  steady  state  is  as  low  as  3  —  4  handovers  per  minute. 


2  Information  Dissemination  in  Satellite  Networks 

2,1  Push-Based  Multi-stage  Information  Dissemination 

Broadcast  data  delivery  is  an  effective  way  to  disseminate  information  to  massive  user  populations. 
The  two  foremost  advantages  of  broadcast  data  delivery  are  scalability  to  large  user  populations, 
and  high  availability  to  geographically  distributed  mobile  users.  However,  users  need  expensive 
and  cumbersome  equipment  to  receive  and  transmit  satellite  signals.  Furthermore,  as  the  amount 
of  information  being  broadcast  increases  average  user  latency  increases  as  well.  Often  users  in 
a  geographical  locality  have  similar  interests,  which  can  be  better  served  by  employing  a  local 
broadcast  schedule.  In  this  context,  a  two  stage  satellite-terrestrial  wireless  broadcast  system 
can  provide  more  efficient  service  in  terms  of  lower  average  user  latency,  and  cheaper  and  more 
convenient  user  equipment.  In  such  a  system,  main  server  broadcasts  information  via  satellite  to  the 
geographically  distributed  local  ground  stations.  Every  ground  station  has  limited  buffer  capacity 
to  store  the  data  broadcast  by  the  satellite.  According  to  their  buffer  content  and  the  interests  of 
their  users  local  stations  deliver  information  via  terrestrial  wireless  channel. 

Main  contribution  of  this  research  initiative  is  the  identification  and  treatment  of  the  issues 
that  arise  in  novel  multi-stage  broadcast  scheduling-caching  applications.  Specifically  in  our  work, 
we  investigate  the  design  of  the  satellite  schedule,  the  cache  management  policies  at  the  ground 
stations  and  the  terrestrial  schedules  that  arise  in  joint  fashion  in  the  context  of  our  problem.  Our 
results  demonstrate  that  such  a  two  stage  system  is  feasible  and  can  provide  more  efficient  data 
delivery  compared  to  a  single  stage  system. 


2.1.1  System  Model 

Under  the  push-based  broadcasting  approach  as  depicted  in  Figure  11,  a  server  continuously  and 
repetitively  broadcasts  data  to  a  user  communitywithout  any  feedback  about  the  users’  needs  due 
to  the  limited  uplink  communication  capability  from  users  to  the  server.  The  data  is  delivered 
in  terms  of  items,  where  an  item  is  a  fixed  sized  information  unit.  The  items  axe  broadcast  in  a 
particular  order  that  is  called  a  schedule.  Time  is  slotted,  and  dmration  of  each  slot  is  equal  to 
the  transmission  time  of  an  item  on  the  broadcast  channel.  When  a  user  needs  a  certain  item,  it 
monitors  the  broadcast  channel  until  the  desired  item  is  detected.  There  is  some  ^latency  for  the 
fulfillment  of  the  needs  of  a  user,  since  the  server  can  only  broadcast  a  single  item  per  broadcast 
channel  at  any  time  instant.  This  latency  depends  on  the  broadcast  schedule,  as  well  as  the  user 
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request  pattern. 


Figure  12  depicts  the  two-stage  satellite-terrestrial  system  that  is  going  to  be  discussed  in  the 
remainder  of  this  paper.  In  this  setup,  the  main  server  transmits  information  items  via  satellite  to 
the  geographically  distributed  local  ground  stations  with  a  transmission  rate  of  Rs  items  per  unit 
time.  The  ground  stations  have  a  limited  local  buffer  capacity  of  C  items.  Each  ground  station 
retransmits  the  items  stored  in  the  local  cache  to  a  local  user  community  with  transmission  rate  Rg 
items  per  unit  time,  where  Rg  =  r  •Rg,  and  r  is  a  positive  integer.  We  investigate  the  problem  of 
minimizing  the  overall  mean  user  latency  by  developing  intelligent  satellite  and  terrestrial  broadcast 
schedule  and  local  server  cache  management  policies.  In  this  context,  the  terrestrial  schedule  should 
try  to  niatch  the  user  access  characteristics  as  much  as  possible,  given  the  satellite  schedule  and 
the  available  items  saved  in  the  cache.  The  cache  management  policy  should  in  turn  store  the 
items  broadcast  by  the  satellite  that  are  expected  to  be  broadcast  by  the  terrestrial  schedule.  The 
satellite  schedule  should  deliver  the  items  in  the  order  that  the  underlying  cache  management  and 
terrestrial  scheduling  policies  can  minimize  the  local  average  user  latency. 


2.1.2  Joint  Cache  Management  and  Terrestrial  Broadcast  Scheduling 

Even  if  we  use  naive  cache  management  and  terrestrial  scheduling  techniques,  the  resulting  perfor¬ 
mance  improvement  with  the  two-stage  systems  is  significant.  Let’s  illustrate  this  with  the  following 
example.  Assume  that  there  are  four  groups  of  users  that  are  interested  in  only  four  items.  Each 
item  is  requested  with  the  same  probability.  Neither  of  the  items  requested  by  the  users  of  differ¬ 
ent  groups  are  the  same.  Therefore,  the  main  server  broadcasts  sixteen  items  to  the  local  ground 
stations  of  these  four  groups.  Assume  that  the  satellite  schedule  is  uniform,  that  is,  these  sixteen 
items  are  broadcast  sequentially  and  periodically.  Each  local  groimd  station  can  store  only  two 
items  at  any  given  time  instant  and  the  satellite  channel  is  twice  as  fast  as  the  terrestrial  channel. 
An  intuitive  and  simple  cache  management  procedure  is  to  “store  any  new  item  of  interest  received 
from  the  satellite  schedule  and  remove  the  item  that  is  the  oldest  in  the  cache.”  Under  this  policy, 
the  cache  content  evolves  as  depicted  in  Figure  13.  An  equally  simple  terrestrial  scheduling  policy 
“broadcasts  at  a  given  instant  the  item  residing  in  the  cache  that  is  broadcast  least  recently.”  Un¬ 
der  this  policy  the  terrestrial  schedule  evolves  as  depicted  in  Figure  13.  It  is  easy  to  see  that  the 
resulting  terrestrial  schedule  has  a  mean  response  time  of  4.25  satellite  slots.  Instead  of  two  stage 
broadcast  data  delivery  if  a  single-stage  broadcast  scheduling  policy  was  deployed,  the  single-stage 
broadcast  schedule  would  have  to  service  all  four  user  groups  and  the  best  schedule  would  then 
be  the  uniform  satellite  schedule  that  is  also  depicted  in  Figure  13,  The  mean  response  time  of 
this  schedule  is  8  satellite  slots.  The  improvement  in  the  average  user  latency  by  employing  the 
two-stage  system  is  almost  50%. 

If  all  items  that  are  requested  by  the  users  in  a  geographical  locality  are  available  at  the  local 
ground  station,  then  there  is  no  need  for  cache  management,  and  the  problem  reduces  to  single  stage 
broadcast  scheduling  problem,  which  is  extensively  studied.  In  these  studies,  it  has  been  noticed 
that  in  any  broadcast  scheduling  system,  there  are  two  quantities  related  to  each  item  i  that  affect 
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the  scheduling  decision  at  slot  n.  These  are,  elapsed  time  Wi{n)  since  the  last  transmission  of 
item  i  and  Ai  the  request  generation  rate  for  item  i.  The  likelihood  of  item  i  being  transmitted 
at  n  increases  with  Ai  and  Wi{n).  In  [10],  Su,  Tassiulas  and  Tsotras  considered  a  set  of  policies 
for  a  single-stage  system,  where  the  broadcast  scheduling  is  determined  based  on  priority  indices 
of  the  items.  The  priority  index  of  item  i  is  the  product  X]wi{n)y  where  7  is  an  exponent  that 
reflects  relative  importance  of  Ai  versus  Wi{Ti).  It  is  shown  that  7  =  0.5  gives  the  best  performance 

among  these  priority  index  policies.  7  =  0.5  refers  to  a  special  case,  where  the  index  is 

Note  that  is  the  aggregate  expected  delay  experienced  by  item  i  requests  since  the  last 

transmission  of  item  i  before  slot  n.  Su  et  al.,  showed  that  this,  so  called  Mean  Aggregate  Delay 
(MAD)y  algorithm  gives  performances  close  to  the  lower  bound. 

The  single  stage  and  the  proposed  two-stage  systems  differ  in  the  availability  of  items  for 
scheduling  at  a  given  slot.  In  two-stage  system,  only  the  items  that  are  stored  in  the  cache  of  the 
local  ground  station  are  available  for  scheduling.  Intuitively,  we  want  the  terrestrial  schedule  of  a 
local  station  in  the  two-stage  system  be  same  as  the  schedule  in  a  single-stage  system  that  is  serving 
exclusively  the  local  user  group.  In  order  to  achieve  this  goal  we  propose  to  use  the  MAD  scheduling 
policy  as  the  local  terrestrial  scheduling  policy  and  develop  good  cache  management  policies  that 
work  on  top  of  the  MAD  scheduling  policy.  A  good  cache  management  policy  should  try  to  feed 
the  scheduler  with  the  items  that  are  going  to  be  broadcast  in  the  coming  slots  according  to  the 
MAD  policy  as  if  all  items  are  available  to  the  scheduler.  We  now  describe  two  cache  management 
algorithms  that  comply  to  this  intuition. 

TWT  Algorithm 

Due  to  the  lack  of  uplink  request  channel,  the  only  information  about  the  user  accesses  that  the 
local  stations  can  determine  at  a  given  slot,  is  the  expected  number  of  users  waiting  for  service. 
The  expected  user  backlog  for  item  i  observed  by  an  arbitrary  local  station  at  satellite  slot  n  is 
defined  as  Xi{n)  =  XiWi{n),  The  evolution  of  the  expected  user  backlog  for  item  i  at  a  local  station 
is  depicted  in  Figure  14.  The  aggregate  average  delay  of  item  i  at  slot  k  is  defined  as  ^Xi{wi{k))^ ^ 
and  it  is  the  striped  area  between  to  and  k  in  Figure  14. 

Assume  that  at  satdlite  slot  k  item  i  is  broadcast  by  the  satellite.  The  local  station  has 
two  options;  to  fetch  item  i  and  replace  another  item  in  the  cache,  or  to  leave  the  cache  content 
unchanged.  If  item  i  is  not  fetched  at  k,  the  earliest  time  it  can  be  scheduled  at  the  terrestrial 
schedule  is  at  satellite  slot  k  +  li{k)^  i.e.  when  item  i  is  re-transmitted  by  the  satellite.  At  that 
instant  the  expected  item  i  backlog  in  terms  of  satellite  slots  will  be  given  by,  bi{k)  =  Xi(k+li(k))  = 
^Xi[wi{k)  +  li{k)]^.  Observe  that  bi{k)  refers  to  the  shaded  area  in  Figure  14.  The  mean  aggregate 
delay  depends  on  the  evolution  of  the  expected  backlog  of  items,  bi{k).  Thus,  higher  is  the  bi{k) 
higher  is  the  mean  aggregate  delay.  TWT  algorithm  fetches  an  item,  if  the  total  waiting  time  of 
that  item  at  its  next  satellite  re-transmission  is  high. 
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ETT  Algorithm 


One  weakness  of  TWT  algorithm  is  that  it  is  a  myopic  algorithm  in  the  sense  that  it  does  not  take 
into  consideration  the  future  desired  terrestrial  schedule.  Unlike  TWT  algorithm,  ETT  algorithm 
first  estimates  the  future  evolution  of  the  terrestrial  schedule  at  an  arbitrary  satellite  slot  k.  Ac¬ 
cording  to  this  estimate,  ETT  algorithm  then  checks  whether  item  i  E  B{k)  is  broadcast  or  not 
at  terrestrial  channel  prior  to  item  i’s  first  re-transmission  at  the  satellite  channel  after  the  fcth 
slot.  If  it  is  not  broadcast,  then  there  is  no  need  to  keep  item  i  in  the  cache.  Else,  eTT  algorithm 
calculates  the  reward  of  keeping  item  i  in  the  cache  at  slot  k.  Let  ti  correspond  to  the  satellite 
slot  the  terrestrial  transmission  of  item  i  is  predicted  to  take  place,  and  let  t2  correspond  to  the 
satellite  slot  the  first  satellite  re-transmission  of  item  i  after  fcth  slot  takes  place  (see  Figure  15). 
If  item  i  is  not  stored  at  fc,  the  earliest  time  item  i  can  be  scheduled  at  the  terrestrial  schedule  is 
at  ^2*  Item  i  is  not  supposed  to  be  broadcast  at  the  terrestrial  channel  until  ti,  so  in  fact  there  is 
no  benefit  of  keeping  item  i  stored  in  the  cache  until  ti.  Thus,  the  reward  of  keeping  item  i  in  the 
cache  is  only  the  additional  waiting  time  from  ti  to  t2  that  is  eliminated  by  the  broadcast  of  item 
i  at  ti.  This  reward  is  illustrated  as  di{k)  in  Figure  15.  Hence,  ETT  algorithm  stores  an  item  i  at 
slot  fc,  if  the  reward  of  keeping  item  i,  di(fc),  is  sufficiently  large. 

In  Figure  16  we  demonstrate  the  performance  of  TWT  and  ETT  with  respect  to  a  lower  bound 
that  we  derived  in  [8].  In  this  example,  the  probability  that  item  i  is  requested  is  given  by  Zipf  1 
distribution,  i.e.  Pr{i)  =  where  c  is  a  normalization  constant  [11].  ETT  outperforms  TWT  with 
quite  a  margin.  This  is  because  ETT  estimates  the  future  terrestrial  schedules  and  tries  to  buflFer 
items  from  satellite  schedule  to  satisfy  this  estimated  schedule.  More  results  can  be  found  in  [8]. 


2.1.3  Satellite  Schedule  Design 


In  general,  the  users  in  different  local  stations  have  different  interests,  i.  e.  different  access  probabil¬ 
ity  distributions.  Therefore  a  satellite  schedule,  that  is  suitable  for  a  user  group,  may  result  in  the 
worst  possible  performance  for  another.  Joint  design  of  the  satellite  schedule,  cache  management 
strategies,  and  terrestrial  schedules  may  provide  a  better  overall  performance.  We  assume  that 
one  of  the  proposed  cache  management  algorithms,  and  the  MAD  scheduling  policy  are  employed 
together  at  the  local  stations.  In  particular,  we  assume  that  the  ETT  algorithm  is  used  as  the 
cache  management  algorithm.  We  consider  off-line  design  of  the  satellite  schedule.  As  a  heuristic, 
we  consider  the  case  where  the  satellite  schedule  is  determined  by  the  MAD  policy.  According 
to  this  problem  setup  the  only  question  that  remains  to  be  answered  is,  which  access  probability 
distribution  should  be  used  to  determine  the  satellite  schedule  that  minimizes  the  overall  average 
response  time  of  the  users.  We  consider  two  cases:  aggregate  probability  distribution  where  all  user 
groups  have  equal  weight,  and  WSAT  where  the  weight  of  each  user  group  is  inversely  proportional 
to  its  local  buffer  size.  Following  table  demonstrates  the  performance  of  these  approaches  on  an 
example. 
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MRT 

Group 

Distribution 

Cache  size 

WSAT 

Aggregate 

1 

zl(l,20,100) 

Cl  =  5 

24.48 

28.31 

2 

zl(20,20,100) 

C2  =  5 

24.83 

28.30 

3 

zl(40,20,100) 

C3  =  5 

25.64 

30.20 

4 

zl(60,20,100) 

C4  =  10 

24.65 

21.72 

5 

zl(80,20,100) 

C5  =  20 

8.07 

8.07 

Overall 

21.534 

23.32 

Table  1:  Example 

2.2  Efficient  Resource  Allocation  for  Pull-Based  Information  Distribution 

As  it  is  the  case  in  the  Internet,  users  often  pull  information  rather  than  waiting  for  the  required 
information  to  be  pushed  to  them.  As  discussed  previously,  many  users  lack  a  high  speed  network 
connection.  Users  usually  have  a  connection  to  a  local  network,  where  a  local  gateway  station  is 
connected  to  the  backbone  network.  In  rural  areas,  satellite  link  appears  to  be  a  feasible  method 
of  connection  for  these  gateways.  However,  satellite  bandwidth  is  expensive  and  it  will  be  cost- 
effective  to  serve  the  requests  locally.  Thus,  local  gateways  usually  employ  caches  to  store  the 
popular  objects.  Until  now,  system  architecture  is  similar  to  the  system  discussed  in  the  previous 
section.  However,  we  further  consider  the  case  where  local  gateways  that  are  in  close  proximity  to 
each  other  are  connected  via  a  terrestrial  link  so  that  the  user  requests  that  cannot  be  satisfied  at 
one  gateway  are  directed  to  another  gateway.  Such  a  flexibility  to  forward  user  requests  among  the 
gateways  allows  user  requests  to  be  satisfied  locally  and  most  probably  with  much  lower  latency. 
The  architecture  is  summarized  in  Figure  17. 

We  consider  a  business  model  in  which  the  primary  customers  of  the  local  gateways  are  the 
publishers  (origin  servers).  The  objective  of  the  publisher  is  to  maximize  its  net  benefit.  The  net 
benefit  is  the  utility  a  publisher  receives  by  delivering  its  content  to  the  users  with  delay  I>,  minus 
the  cost  of  the  resources  required  to  ensure  that  delay.  A  cost  arises,  because  the  gateway  has 
limited  resources  such  as  bandwidth,  cache,  etc.  and  uses  prices  to  induce  efficient  allocation.  The 
objective  of  a  gateway  is  to  maximize  its  revenue,  while  satisfying  the  requested  QoS.  Publishers 
act  non-cooperatively  and  try  to  maximize  their  own  benefit  regardless  of  the  rest  of  the  servers. 

Content  Delivery  problem  can  be  separated  into  two:  content  distribution  and  user  request 
routing  sub-problems.  Given  user  request  arrival  rates  to  the  surrogates,  the  distribution  sub¬ 
problem  determines  the  best  dissemination  strategy,  i.e.  the  set  of  objects  disseminated  to  each 
surrogate,  maximizing  its  net  benefit.  Meanwhile,  given  the  surrogate  cache  content,  the  request 
routing  sub-problem  determines  the  best  user  request-routing  strategy,  i.e.  the  set  of  user  requests 
served  by  each  surrogate,  again  maximizing  the  net  benefit  of  the  publishers.  Although  these 
problems  are  jointly  related,  for  the  purpose  of  designing  a  practical  algorithm,  we  envision  an 
iterative  scheme,  in  which  each  sub-problem  is  solved  separately.  In  this  scheme,  the  distribution 
sub-problem  determines  the  dissemination  strategy  according  to  the  arrival  rates  as  determined  by 
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the  request-routing  sub-problem.  The  request  routing  sub-problem  determines  the  routing  strategy 
according  to  the  cache  content  as  determined  by  the  distribution  sub-problem.  These  sub-problems 
update  their  decisions  iteratively  according  to  the  output  of  the  other.  We  expect  that  the  solution 
to  the  overall  content  delivery  problem  is  (near-)  optimal,  if  the  solutions  to  these  two  sub-problems 
are  (near-)  optimal.  The  formulation  and  the  aim  of  these  two  problems  closely  resemble  each  other. 
In  the  following,  we  describe  the  results  for  the  distribution  sub-problem  with  the  understanding 
that  similar  results  apply  to  the  routing  sub-problem  as  well.  For  fmrther  details  on  the  routing 
sub-problem  we  refer  the  readers  to  [12], 


2.2.1  System  Model 


The  objective  of  the  distribution  sub-problem  is  to  maximize  the  reduction  in  average  user  latency 
by  the  use  of  surrogate  caches,  while  the  total  cost  of  service  is  less  than  the  highest  investments  the 
publishers  are  willing  to  make  for  the  caching  resources.  Figure  18  illustrates  the  network  set-up 
that  we  are  interested  in  this  section.  There  are  several  LANs  from  which  the  user  requests  originate. 
Every  user  is  interested  in  one  or  more  of  the  objects.  The  user  requests  are  first  intercepted  by 
the  request-routing  sub-system.  The  request-routing  subsystem  then  forwards  the  user  request  to 
an  appropriate  surrogate,  perhaps  the  one  with  the  highest  probability  of  delivering  the  requested 
object  with  minimum  latency.  For  the  distribution  sub-problem,  we  assume  the  routing  decision  is 
taken,  and  the  user  request  arrival  rates  to  each  surrogate  are  determined.  Surrogates  are  located 
between  the  user  networks  and  the  origin  servers.  Thus,  users  are  always  at  most  two  hops^  away 
from  the  content.  A  user  request  is  first  checked  at  the  corresponding  surrogate.  If  the  requested 
object  is  available  at  the  surrogate,  request  is  immediately  served.  Otherwise,  request  is  forwarded 
to  the  web  server,  where  the  object  originally  resides. 

Surrogates  have  limited  cache  size,  and  the  cache  is  shared  among  the  publishers.  Surrogates 
charge  the  publishers  for  the  portion  of  the  cache  their  content  occupies.  We  assume  that  the 
surrogates  does  not  cooperate  -  other  than  announcing  their  resource  prices-  probably  due  to 
high  messaging  and  processing  overheads  associated  with  cooperative  caching.  Instead,  they  act 
non-cooperatively  with  the  objective  of  maximizing  their  individual  revenues.  Surrogates  compete 
among  each  other  to  store  the  publishers  content.  Publishers  do  not  collaborate  either  and  they 
try  to  purchase  as  much  cache  space  as  possible  with  the  objective  of  maximizing  their  net  benefit. 
Assume  that  there  are  I  different  publishers  (origin  servers)  and  J  surrogates  present  in  the  network. 

User  requests  arrive  from  N  different  user  LANs.  Let  denote  the  total  request  arrival  rate  from 

LAN  n  at  surrogate  j  for  the  content  in  the  ith  publisher.  Let  Aj  =  A^’”  be  the  total  arrival  rate 

to  the  surrogate  j  for  the  content  in  publisher  i.  The  user  interest  in  the  objects  of  the  publishers 
is  distributed  according  to  Zipf  distribution  [11].  ' 

Let  be  the  investment  of  the  ith.  publisher  in  the  jth  surrogate.  Let  Bi  =  Y2j  be  the 
^Notice  that  user  packets  may  traverse  multiple  network  elements/links  to  reach  the  server. 
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total  investment  of  the  ith  publisher.  It  is  assumed  that  the  information  stored  in  the  servers  is 
continuous  and  can  be  replicated  continuously  to  a  surrogate.  Total  information  available  at  the 
piiblisher  i  is  Xz*  Publisher  replicates  its  most  popular  part  of  the  content  to  the  surrogates  so  that 
the  cache  hit  probability  is  maximized.  Assuming  that  Ci  units  of  cache  space  is  allocated  to  the 
publisher,  the  probability  that  an  incoming  user  request  is  satisfied  at  the  surrogate  is  given  by, 

Pr(hitlQ)  =  *  q[x)  dx  =  ^  • 

Let  pj  denote  the  price  of  the  unit  cache  space  in  surrogate  j.  Let  the  pricing  policy,  p  = 
(Pi,P2j  •  •  •  jPj)>  denote  the  set  of  unit  cache  space  prices  of  all  the  surrogates  in  the  network.  Let 
dij  denote  the  additional  average  delay  a  user  request  forwarded  from  surrogate  j  to  the  origin 
server  of  publisher  i  will  experience.  Note  that  if  we  consider  geosynchronous  earth  orbit  satellites, 
we  can  assume  that  dij  is  approximately  the  same  for  all  i,  j.  However,  in  this  study  we  consider 
a  more  general  system,  where  a  low-earth  orbit  satellite  network  (LEO)  may  be  used  to  connect 
the  origin  servers  and  the  gateways.  LEO  satellites  have  propagation  delays  comparable  to  the 
terrestrial  networks  and  thus  the  locations  of  the  origin  servers  and  the  gateways  have  effect  on  the 
observed  delays. 


2.2.2  Optimal  Publisher  and  Surrogate  Strategies 

We  first  consider  the  optimal  strategy  for  the  publishers.  Let  be  the  cache  space  allocated  to 
publisher  i  in  surrogate  j.  Hi  th  publisher’s  investment  in  the  j  th  surrogate  is  then  the  total 

cache  space  allocated  to  the  content  of  publisher  i  in  surrogate  j  is  =  ■^.  The  average  reduction 
in  the  user  delay  or  equivalently  the  average  net  benefit  that  publisher  i  generates  by  investment 

j  .  Define  as  the  gain  factor  for  publisher 

i  from  surrogate  j.  The  utility  function,  i.e.  the  total  additional  average  benefit,  Ui{xi),  of  publisher 
i  is  Ui{xi)  —  Ylj=iPi  .  For  a  given  pricing  policy  p  the  publisher  optimization  problem 

(S)  can  be  given  as: 

(S)  Ui{xi) 

J 

subject  to  ^2  ^iPj  ~ 
j=i 

Since  Ui{xi)  is  a  concave  function,  and  the  constraint  set  is  compact,  there  exists  a  unique 
solution  to  (S). 
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Lemma  1 


£i 

P} 


Efc=l  Pk 


Pk 


is  the  unique  optimal  solution  to  the  optimization  problem  (S). 


We  now  consider  the  optimal  pricing  strategies  of  the  surrogates  maximizing  their  revenues. 
Let  =  (pi,P2j  •  •  •  •  •  -  yPj)  be  the  set  of  unit  cache  space  prices  of  all  the  surrogates 

in  the  network  except  the  jfth  one.  We  assume  that  there  is  no  collaboration  among  the  surrogates, 
and  teach  surrogate  tries  to  maximize  its  revenue  non-cooperatively. 

Lemma  2  Surrogate  revenue  rj{pj) 

when  pj  satisfies  Yli=ix\{pj) 
surrogate  cache  space. 


=  ^  maximized  under  a  given  fixed  pricing  policy 

=  Cjj  i.e.  when  the  total  publisher  demand  is  equal  to  the 


Until  now,  we  discussed  the  optimal  strategies  of  the  publishers  and  the  surrogates  given  that 
system  is  at  a  steady  state.  However,  we  have  not  discussed  whether  such  a  steady  state  exists. 
Notice  that  when  a  surrogate  re-evaluates  its  pricing  policy  according  to  the  pricing  policies  of  the 
rival  surrogates,  the  remaining  surrogates  will  do  the  same.  For  each  different  pricing  policy  the 
publishers’  optimal  investments  will  be  different  as  well. 

In  order  to  understand  the  behavior  of  the  surrogates,  we  model  the  two-stage  publisher- 
surrogate  system  as  a  non-cooperative  game.  In  this  publisher-surrogate  distribution  game^  r(  J,  »S,  P), 
the  players,  J,  are  the  surrogates,  the  strategy  set  Sj  for  a  surrogate  j  is  given  by  the  surrogate’s 
unit  cache  space  price  and  the  payoff  function  Pj{s)  of  each  surrogate  j  is  given  by  the  profit  of  the 
jth.  surrogate.  This  system  is  similar  to  the  Cournot  oligopoly  discussed  in  the  economics  literature. 
Assume  that  each  surrogate  has  a  fixed  cost  for  its  cache,  but  has  no  control  over  the  size  of  the 
cache,  i.e.  the  size  of  the  cache  is  determined  prior  to  the  system  operation.  We  first  show  that  this 
game  has  a  Nash  equilibrimn  solution,  where  no  sxmrogate  has  incentive  to  change  unilaterally  its 
strategy,  since  each  surrogate  maximizes  its  own  individual  payoff  given  the  strategies  of  others. 


Theorem  1  The  non-cooperative  publisher-surrogate  distribution  game  r(  J,  5,  P)  has  at  least  one 
Nash  Equilibrium  solution. 


Pareto  optimality  is  the  relevant  criteria  in  a  multi-objective  problem  setting.  At  the  pareto 
optimum  solution,  one  can  find  no  other  feasible  solution  that  increases  some  objectives  while  not 
decreasing  at  least  another  objective.We  have  shown  in  [9]  that  if  the  equilibrium  is  unique,  then  the 
outcome  of  the  non-cooperative  game  is  the  optimal  solution  to  the  individual  revenue  maximization 
problems  of  the  publishers.  If  there  are  multiple  equilibria,  however;  the  resulting  cache  allocations, 
X,  are  only  locally  optimum.  Unfortunately,  there  are  often  multiple  Nash  equilibria  and  depending 
on  the  initial  prices  as  well  as  price  update  strategies  the  outcome  of  the  game  may  not  always  be 
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the  optimal  solution.  In  [9],  we  discuss  a  special  case  of  the  surrogate  cache  allocation  problem, 
where  the  delay  between  each  publisher-surrogate  pair  is  the  same  and  user  request  arrival  rates  to 
each  surrogate  and  Zipf  distributions  for  each  publisher  are  identical.  For  this  case,  we  determine 
the  condition  for  which  unique  equilibrium  exists. 


2.2.3  Selection  of  Publisher  Investments 


Until  now,  we  assumed  that  publisher  investments  are  fixed  and  known.  We  now  discuss  several 
possibilities  for  determining  the  investment  amounts.  Consider  the  following  optimization  problem 
(P).  The  objective  of  (P)  is  to  maximize  the  total  publisher  utility  given  the  limited  cache  sizes. 
We  define  the  solution  of  (P)  as  the  system  optimum  solution. 

(P)  max  y]l7i(xi) 

Kki  i 

subject  to  ^  xl  <Cj,  j  =  1,2, . . .  ,J. 
i 

This  problem  can  be  solved  by  considering  the  dual  problem  and  using  the  gradient  projection 
method.  One  important  feature  of  this  method  is  the  fact  that  it  can  be  implemented  in  a  dis¬ 
tributed  fashion.  The  solution  of  this  problem  dictates  the  publishers  to  pay  the  cost  associated 
with  the  system  optimum  solution.  Notice  that  a  server  may  be  individually  better  off  by  investing 
less  than,  the  amount  dictated  by  the  system  optimum  solution.Thus,  a  system  optimum  solution 
can  only  be  achieved  by  the  distributed  publisher-surrogate  game  by  enforcing  the  publishers  to 
pay  the  charges  associated  with  the  system  optimum  solution.  Without  such  enforcement,  we  en¬ 
counter  a  different  optimization  problem  where  the  investment  amount  is  also  a  design  parameter. 
In  the  resulting  non-cooperative  game,  publishers  select  their  investment  amounts  to  maximize 
their  profits. 

Consider  a  publisher-surrogate  game  where  a  publisher  selects  an  investment  maximizing  the 
“net”  publisher  benefit.  That  is,  the  difference  of  the  additional  average  benefit  generated  from 
using  surrogates  and  the  cost  of  using  them  is  maximized.  Let  T  be  the  duration  of  time  a  publisher 
can  rent  the  caching  space  from  a  surrogate.  Then  for  a  given  set  of  surrogate  prices  p  the  publisher 
optimization  problem  can  be  given  as, 


max  Ui{xi) - —  (7) 

subject  to  ^2  oc^pj  <  Bi- 

j 


We  developed  a  suboptimal  investment  strategy  for  the  publishers  by  assuming  that  the  change 
in  the  equilibrium  prices  is  small  when  the  change  in  the  investments  is  small  as  well.  For  a  given 
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investment  Bi,  the  optimal  cache  allocation  is  given  by  Xij{Bi)  —  Assume  that  pj  are 

the  set  of  prices  at  the  equilibrium.  As  for  the  key  assumption  of  our  derivation,  we  assume  that 
for  a  small  change  in  investment  amounts,  the  change  in  the  equilibrium  prices  is  small  as  well. 
Specifically,  in  the  vicinity  of  the  set  of  investment  amounts  Y^kPfk/pk  can  be  considered  as 
constant  This  assumption  is  verified  through  numerical  studies.  Assuming  that  this  condition  is 
satisfied,  the  optimal  investment  amount  maximizing  the  net  benefit  is  calculated. 


2.2.4  Numerical  Results 

Figures  19  and  20  demonstrate  the  performances  of  the  suggested  approaches.  Let  C/  = 

Prom  the  definition  of  Cf  and  assuming  that  Wi{-)  is  the  same  for  all  publishers,  one  can  notice  that 
by  varying  Cf ,  we  basically  change  the  delay  between  surrogate  j  and  publisher  i.  As  illustrated  in 
Figure  19,  as  the  skewness  of  the  system  increases  the  performance  of  game-theoretic  algorithm  gets 
better  compared  to  the  conventional  algorithm.  However,  this  improvement  is  larger  for  low  values 

of  C/ •  This  means  that  surrogates  need  not  be  very  close  to  the  users.  We  may  place  surrogates  at 
locations  at  somewhat  greater  distance  and  thus  service  more  users  with  less  number  of  surrogates 
with  sufficiently  good  system  performance. 

In  Figure  19  the  a  priori  fixed  investments  are  81=82  =  10.  The  system  optimum  solution 
maximizes  the  total  publisher  utilities.  As  expected  we  observe  that  the  total  utility  achieved 
by  the  investments  associated  with  the  system  optimum  solution  is  indeed  higher  than  the  fixed 
investments.  It  is  interesting  to  note  however,  while  the  system  optimum  investment  maximizes 
the  total  of  publisher  utilities,  it  reduces  some  of  the  publishers  individual  utilities.  Furthermore, 
the  investment  amounts  leading  to  the  system  optimum  are  usually  very  high,  reducing  any  benefit 
received  from  the  use  of  the  caches  in  the  first  place.  A  short-sided  objective  of  maximizing  total 
system  utility  is  thus  seems  not  implementable,  when  the  system  is  formed  by  selfish  non-cooperative 
agents. 

.  In  Figure  20,  we  investigate  the  performance  of  our  sub-optimal  investment  algorithm  that 
prevents  the  abovementioned  problem.  We  compare  our  results  with  the  game  theoretical  algorithm 
with  fixed  investments,  where  the  investments  are  10  units.  The  system  setup  is  similar  to  the 
previous  numerical  examples,  and  the  improvement  is  considered  for  varying  gain  factor.  We  also 
depicted  the  reduction  in  the  total  publisher  utility  as  compared  to  the  optimal  system  utility.  The 
cache  rental  period  is  100  time  units.  In  this  example,  we  observe  that  the  improvement  in  the 
publisher  profits  is  significant  when  the  publishers  are  allowed  to  optimally  vary  their  investments. 
Furthermore,  the  decrease  in  the  system  efficiency  is  small  as  compared  to  the  increase  in  the  profit. 


21 


2.2.5  An  Optimization  Based  Resource  Allocation  Algorithm 


1  t 


■f 


The  results  given  in  the  previous  sections  suggest  that  we  may  use  a  price-directed  market-based 
distributed  algorithm  for  solving  the  two-stage  publisher-surrogate  cache  resource  allocation  prob¬ 
lem.  We  consider  the  following  algorithm  for  this  purpose:  An  initial  set  of  prices  is  announced 
to  the  publishers.  The  publishers  determine  their  resource  (cache)  demands  according  to  these 
prices  as  well  as  the  request  rates,  and  the  observed  delays  from  the  surrogates.  The  publishers 
request  these  resources  from  the  surrogates.  Prices  are  then  iteratively  changed  to  accommodate 
the  demands  for  resources  until  the  total  demand  equals  to  the  total  amount  of  resource  available. 

In  [9],  we  show  that  this  algorithm  converges  to  the  Nash  equilibrium  solution. 
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(a)  Carrier  with  one  traffic  burst  offset 
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n?)  Carrier  with  several  traffic  burst  offsets 
Figiure  1:  Comparison  of  the  cases  of  single  and  multiple  burst  offset  values. 


Satellite 


Figure  2:  Different  delay  classes  within  a  satellite  footprint. 
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Figure  3:  Relative  position  of  delay  classes  and  beams  in  the  satellite  footprint. 
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Figure  4:  Illustration  of  first  step  of  delay  class  determination  algorithm. 
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Figure  5:  A  beam  with  two  delay  classes  and  the  corresponding  overlap  region. 


Figure  6:  Schematic  representation  of  Path  Selection  procedure. 
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Figure  7:  Demonstration  of  beam  handover  from  the  beam  located  in  the  satellite  nadir  (C)  to  a 
neighboring  beam  (F). 


BLOCKING  PROBABILITY  FOR  AN  EDGE  BEAM  Wrm  12  CARRIERS 


BLOCKING  PROBABIUTY  FOR  AN  EDGE  BEAM  WITH  4S  CARRIERS 


Figure  8:  Blocking  rates  with  and  without  the  carrier  grouping  method,  for  an  edge  beam  with  12 
and  '45  carriers. 
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AGGREGATE  SATELUTE.  BEAM  AND  DELAV  CLASS  HANDOVER  RATES  DURING  BUSY  HOUR  AT  (5.17)  ENTIRE  GLOBE  :  SATELLITE,  BEAM  AND  DELAY  CLASS  HANDOVER  RATES  DURING  BUSY  HOUR 


Figure  9:  Satellite,  beam  and  delay  class  handover  rates:  Handover  rates  at  a  square  region  of  15° 
longitude  and  latitude  distance  from  the  GW,  with  A  =  4.87  calls/sec  and  handover  rates  for  all 
calls  served  by  a  GW. 


Figure  10:  Satellite  and  beam  handover  rates  under  UT  position  or  maximum  beam  residence  time 
handover  criterion. 
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Figure  ll;  A  Single-stage  Broadcast  Data  Delivery  System  in  a  Wireless  Conai 
ment 


Main  server  Local  Ground  station  User  Community 

Figure  12:  Multistage  Information  Distribution 


Figure  13:  The  cache  management,  and  terrestrial  schedule  at  one  of  the  four  ^ 
users  of  this  ground  station  are  equally  interested  in  items  1,  5,  9,  13.  The  sate 
as  fast  as  the  terrestrial  channel. 
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Figure  14:  The  evolution  of  the  expected  user  backlog  of  item  i  observed  by  a  local  station.  The 
upper  side  of  the  u-axis  refers  to  the  satellite  transmissions,  while  the  lower  side  refers  to  the 
terrestrial  transmissions. 


First  estimated  terrestrial 
broadcast  of  item  i 


Figure  15:  Illustration  for  ETT  algorithm.  The  upper  side  of  the  u-axis  refers  to  the  satellite  trans¬ 
missions  while  the  lower  side  refers  to  the  terrestrial  transmissions.  ti  denotes  the  first  estimated 
terrestrial  transmission  slot  of  item  i  after  kth.  slot. 
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Figure  16:  The  performances  of  TWT  and  ETT.  The  satellite  schedule  is  uniform  with  100  items 
to  broadcast,  and  the  user  access  probability  distribution  is  Zipf  1  distribution  over  100  items 
broadcast. 


LfOcal  User  Networks 


Figure  17:  A  general  multi-stage  satellite  information  retrieval  system. 
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Figure  18:  Content  delivery  system. 
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Figure  19:  The  percentage  of  improvement  of  using  the  game-theoretic  algorithm  with  varying 
investments  over  the  conventional  caching  methods.  ai  =  0.1,  a2  =  0.4.  == 

1,  V(i,i)  ^  (1, 1),  (2,3).  A}  =  =  A?  =  A3  =  3.  X{  =  l,V(i,  j)  ^  (1, 1),  (2, 1),  (1,3),  (2,3). 


Figure  20:  The  percentage  of  improvement  of  usin  g  the  game-theoretic  algorithm  with  varying 
investments  over  the  conventional  caching  methods.  The  comparison  is  made  with  fixed  investments 


where  =  B2  =  10.  ai  =  0.5,  aa  =  0.5.  Cl  =  C|  =  C-  C{  =  l,V(i,y)  ^  (1,1),  (2, 3).  Aj  =  A^  = 
Af  =  Ai  =  3.  Af  =  l,V(i,j)^(l,l),(2,l),(l,3),(2,3). 
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